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Abstract. Calculations of the spectral cloud reflectance over bare soil 
and vegetation are performed in the framework of the asymptotic radia- 
tive transfer theory. A simple equation for the threshold value (THV) of 
the albedo of the underlying surface is derived. The cases of underlying 
surfaces with the albedos smaller than the THV can be treated as if the 
albedo of the surface is equal to zero in satellite cloud remote sensing 
problems.. It was found that the THV increases with the cloud optical 
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1 Introduction 


Spectral reflection and transmission characteristics of turbid 
media are of importance for a number of applications (e.g., 
oceanic and atmospheric optics). Reflection and transmit- 
tance of light by random media can be studied in the frame- 
work of radiative transfer theory, which is based on the 
integro-differential radiative transfer equation.'” Theoreti- 
cal results can be obtained both for an underlying black 
surface (e.g., clouds over ocean in the IR) and for bright 
surfaces (e.g., sand and snow in the visible). 

One can derive the following simple relationship be- 
tween the cloud reflectance R, for the case of black ground 
and the cloud reflectance R for the case of a Lambertian 
ground with albedo A? 


Atty 


R=R,+ ` 
1-Ar 


(1) 


where R,=R(A=0), r is the spherical albedo of a cloud 
over the black surface, and tj is the cloud’s total transmit- 
tance in the direction specified by the solar (j=1) or obser- 
vation (j=2) zenith angle (at A=0). Zenith angles are mea- 
sured from the normal to a scattering layer. 

The measured value of R is often used for the estimation 
of the optical thickness of the medium. In particular, satel- 
lite remote sensing of clouds is performed by fitting mea- 
sured spectral cloud reflectances R(\) (where N is the 
wavelength) to synthetic ones.” The problem is that in the 
case of thin clouds the underlying surface may considerably 
change the registered spectra. Therefore, there is a need to 
take account of the underlying surface reflectance A in 
cloud retrieval algorithms. The value of A is not known a 
priori, which constitutes a major problem for satellite thin- 
cloud retrieval algorithms. However, in some cases (thick 
clouds over dark surfaces, e.g., water) there is no necessity 
to take account of the surface contribution, because it is 


very weak. The corresponding measurements of A(A) can 
be found, e.g., in the work of Krinovy.* 

The main task of this paper is to calculate the boundary 
value of the underlying surface reflectance A“ having the 
property that the reflection function R of the cloud— 
underlying-surface system differs by not more than a given 
amount s from the value of R at A=0 for all ASA’. This 
information is of importance for understanding if the sur- 
face contribution must be taken into account in the retrieval 
algorithm for a given cloud field. The value of s can be 
taken equal to the calibration error of the optical instrument 
used for the measurements. 

We also considered the influence of underlying surface 
on the cloud transmittance. 


2 Reflectance 


The boundary value of A” can be easily found from Eq. (1). 
Namely, it follows that 


A tt 
—— 2 e, (2) 
(1 -A r)R, 


where s=(R-R,)/R,. We obtain from Eq. (2) 


yen (3) 
erR, + tita 

This equation enables the estimation of the value of A” for 
given values of ¢;, r, and R,. One can expect the influence 
of the surface to be less important for optically thick clouds 
[t;—0, r— 1, A*— 1; see Eq. (3)]. In this case the value of 
R, for nonabsorbing media can be written approximately 
as 


Rp = Rs- tKo()Ko(), (4) 
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Fig. 1 Dependence of A’ given by Eq. (3) on cloud optical thickness 
7 for ice (broken line) and water (solid line) clouds for nadir obser- 
vation and a solar zenith angle of 60 deg. 


1 
f= 3 
1.072 + 0.757(1 — g) 


(5) 


Tis the cloud optical thickness, g is the asymmetry param- 
eter, € is the cosine of the solar zenith angle, 77 is the cosine 
of the observation zenith angle, R., is the reflection function 
of a semi-infinite cloud, and Ko(€) is the escape function, 
which can be written as 


3 
Ko(é) = zl + 26) (6) 


for 20.2. The accuracy of Eq. (4) is studied by Kokha- 
novsky and Rozanov,’ who show that Eq. (4) is applicable 
for T=5. It is easy to demonstrate that” r=1—-t, t;=tK(&), 
and t,=tKo(7) for nonabsorbing optically thick turbid me- 
dia. 

The results of calculations of the maximal albedo of the 
underlying surface [according to Eq. (3)] for nadir observa- 
tion and a solar zenith angle of 60 deg for optically thick 
nonabsorbing clouds are shown in Fig. 1 for water clouds 
(solid line, g=0.85) and for ice clouds (dashed line, g 
=0.75) at e=0.05. The value of R, was taken from exact 
radiative transfer calculations for water (R,=0.87) and ice 
(R,=0.92) clouds. It follows that remote sensing of clouds 
with 7=5 over water (A <0.05) can be performed neglect- 
ing the contribution of the surface reflection. Clearly, this is 
not the case for bright surfaces. For instance, the surface 
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Fig. 2 Spectral albedo of vegetation (solid line) and bare soil 
(dashed line) (W. von Hoyningen-Huene, private communication). 


reflectance of snow is typically 0.8 in the visible. For water 
clouds with the optical thickness 70 and below, the contri- 
bution of the snow surface must be taken into account. The 
correspondent optical thickness is equal to 40 for ice clouds 
over snow, due to the smaller transmittance of ice clouds 
than that of water clouds having the same optical thickness. 

Land surface reflectance (outside of complex urban ar- 
eas) can be influenced by a mixture of the spectra of veg- 
etation and bare soil. Some examples are given in Fig. 2. It 
follows from Fig. 2 that the reflectance can reach 0.3 in the 
near IR for the vegetation. This means that cloud retrievals 
over land cannot ignore the contribution from the second 
term in Eq. (1) (see also Fig. 1), especially for thin clouds. 
The results of calculations of the spectral cloud reflectance 
over land using Eq. (4) are shown in Figs. 3 and 4. The step 
around 0.7 um is clearly seen in Fig. 3 at t<10. This 
means that the effect of the surface reflectance is important 
in this case. The surface effects can be neglected at 7> 30, 
which is in agreement with Fig. 1. The results for bare soil 
are shown in Fig. 4. Here we see a slight increase of the 
cloud reflectance with the wavelength, due to the influence 
of the underlying surface, at rs 20. The slope of the curves 
decreases with increasing cloud optical thickness. 

The reflection function results for water and ice clouds 
are shown in Fig. 5 at 7=5. The step in the spectral reflec- 
tance curve R(A) is more pronounced for water clouds be- 
cause water clouds are more transparent (larger values of g) 
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Fig. 3 Reflection function of a water cloud for nadir observation and 
a solar zenith angle equal to 60 deg, at cloud thicknesses equal to 
5, 10, 20, and 30. The spectral albedo of vegetation as shown in Fig. 
2 was used in the calculations. 


to light coming from the ground. So the contribution of the 
surface for water clouds is larger than for ice clouds in the 
visible. 


3 Transmittance 


A similar procedure can be followed to study the influence 
of the ground reflectance on the transmission function T of 
a cloud. This function is given as? 


T=T,+ 5, (7) 


where T,=T(A=0), r, is the cloud plane albedo, and the 
other parameters are the same as in Eq. (1). The value of A” 
for the transmitted light can be obtained from Eq. (3) by 


1.0 
0.8 
0.7 


0.6 


0.5 


0.4 


reflection function 


0.3 


0.2 


0.1 


0.0 
0.3 0.4 05 0.6 07 0.8 0.9 1.0 1.1 


wavelength, nm 


Fig. 4 The same as in Fig. 3, but over bare soil. 
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Fig. 5 Reflection function of water and ice clouds over vegetation 
for nadir observation and a solar zenith angle equal to 60 deg, at 
cloud thickness 5. The spectral albedo of vegetation as shown in 
Fig. 2 was used in the calculations. 


replacing R, with T, and also replacing t with r,. It fol- 

lows that 

Ie 
erT, + trp 


(8) 


Then using the results of the asymptotic theory for optically 
thick nonabsorbing clouds,“ 


T=tK(4)Ko(n), tı =tKo(€), rp=1- tKo(9), (9) 


we can easily derive the following expression for the 
threshold value of A“: 


~ e(1—) 1+ K) 


A (10) 


This gives a maximal ground albedo of the surface, 
which must be used in the calculations of T, with accuracy 
given by the threshold value £. The dependence of A“ on T 
is given in Fig. 6 for ice and water clouds. It is seen that A” 
is smaller for ice clouds. This is easily understood. Indeed, 
ice clouds have smaller values of the asymmetry parameter, 
and therefore ice clouds having the same optical thickness 
reflect more and have a much stronger interaction with the 
ground surface for nonabsorbing channels than do water 
clouds (in the visible, where light absorption processes can 
be neglected). So the ground surface plays a larger role in 
this interaction process. This makes A“ smaller than for 
water clouds (see Fig. 6). 

The dependence of the transmission function on the 
wavelength for underlying vegetation and for bare soil is 
shown in Figs. 7 and 8, respectively. The surface reflec- 
tance spectra are assumed to be the same as shown in Fig. 
2. 

It follows that there is some influence of the underlying 
surface albedo on the transmitted light even for very thick 
clouds (see Fig. 7). In particular, we see that the transmis- 
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Fig. 6 Dependence of A’ [see Eq. (8)] on cloud optical thickness for 
ice (broken line) and water (solid line) clouds at the zenith observa- 
tion (€=1). 


sion for larger wavelengths is somewhat larger due to the 
surface reflectance effects. This will lead to a difference 
between the spectral compositions of the solar light below 
and above the cloud. 
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Fig. 7 Transmission function of water cloud for zenith observation 
(€=1) and a solar zenith angle equal to 60 deg at cloud thicknesses 
equal to 5, 10, 20, 30, and 100. The spectral albedo of vegetation as 
shown in Fig. 2 was used in calculations. 
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Fig. 8 The same as in Fig. 7, but over bare soil. 


4 Conclusions 


We have found that the reflection function of water clouds 
with optical thicknesses larger than approximately 30 can 
be calculated assuming a black underlying surface. This 
peculiarity can be used for the simplification of cloud re- 
trieval algorithms over land in the visible and near IR. In 
particular, it is of advantage to use smaller wavelengths 
(e.g., in the UV) for cloud remote sensing. Then the influ- 
ence of the surface is minimized. 

For smaller values of 7, one must take account of the 
surface reflectance both in forward and in inverse model- 
ing. The situation is much alleviated over water, whose 
light reflectance is generally below 5%. Then the reflec- 
tance of light from clouds with T= 5 can be treated without 
consideration of the surface reflectance. This simplifies re- 
trieval algorithms considerably for extended cloud fields. 
Clearly, the surface contribution cannot be neglected for 
broken-cloud conditions even if the surface albedo is equal 
to zero and the optical instrument observes a ground scene 
with some portions not completely covered by a cloud. 
Then measurements of R(A) are governed by the cloud 
fraction. This complex situation, however, is outside the 
scope of this paper, which deals with the idealized case of 
horizontally and vertically homogeneous clouds. 

Also, the case of diffuse transmitted light has been con- 
sidered. It was found that the ground surface influences the 
transmission spectra even for relatively thick clouds. This 
is in contrast with the reflectance spectrum of an extended 
cloud field. 

Our results suggest that the spectral top-of-atmosphere 
reflectance measured by a spectrometer on a satellite orbit- 
ing a planet can be used to identify the type of the under- 
lying surface underneath (e.g., the presence of vegetation; 
see Fig. 3) even for the case of a cloudy atmosphere (so 
long as the ground albedo is larger than the threshold value 
A” introduced in Sec. 1). 

To conclude, we underline that the theoretical results 
obtained are of importance for several other applications of 
optical engineering, including the paint, paper, and textile 
industries. 
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